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ABSTRACT 


This  paper  Is  concerned  with  the  "Square  Root  Damage 
Law" — an  analytical  method  of  computing  blast  damage  on  target 
complexes.  The  Law  assumes  that  target  value  Is  a continuous 
function  of  position,  and  replaces  a set  of  Individual  weapons 
attacking  a target  by  a weapon  density  function  that  Is  optimized 
to  give  maximum  damage.  This  basic  approach  Is  applied  to 
several  types  of  target  damage  functions  and  methods  of  weapon 
targeting. 

Attention  is  restricted  to  cases  where  the  results  can  be 
presented  in  relatively  simple  analytical  form.  The  derivations 
are  presented  in  some  detail  to  illustrate  tne  nature  of  the 
results  obtained.  Charts  of  darnage  are  presented  for  sets  of 
dimensionless  parameters  that  govern  the  damage  sensitivity  to 
weapon  usage.  As  a practical  application,  they  can  form  analyt- 
ical bases  for  explaining  the  results  obtained  from  detailed 
numerical  calculations  of  the  damage  for  a particular  target 
from  a specific  set  of  weapon  aimpoints. 
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Chapter  I 
INTRODUCTION 


In  assessing  the  effects  of  nuclear  attack  upon  target 
complexes,  attention  often  is  restricted  to  blast  effects,  or 
more  precisely  to  effects  where  the  probability  of  damage  is 
taken  only  as  a function  of  the  distance  from  a weapon,  inde- 
pendent of  direction.  In  many  cases  the  target  complex  is 
assumed  to  consist  of  a set  of  discrete  target  elements;  each 
element  has  an  overall  probability  of  damage  obtained  by  com- 
bining the  probabilities  of  damage  from  separate  individually 
located  weapons;  and  overa.ll  damage  is  calculated  by  summing 
the  effects  on  each  target  element.  While  this  type  of  calcula- 
tion is  the  natural  way  of  assessing  blast  damage  from  a 
particular  attack,  the  results  are  strictly  numerical  and  no 
analytical  insight  is  gained  from  such  a calculation. 

Analagous  to  the  dichotomy  between  particulate  and  continuum 
mechanics,  an  approach  to  casualty  estimation  can  be  developed 
where  target  elements  are  replaced  by  a target  value  as  a con- 
tinuous function  of  position,  and  weapon  locations  are  replaced 
by  a weapon  density  function.  When  the  target  value  function 
J,s  a circular  normal  function,  the  local  fraction  of  damage 
is  an  exponential  function  of  weapon  density,  and  weapons  are 
optimally  targeted  against  value,  a simple  expression  can  be 
obtained  for  overall  damage  to  the  target  complex  as  a function 
of  weapon  usage.  This  expression  has  become  known  as  the 
’’Square  Root  Damage  Law.”  A derivation  of  this  law,  using  a 
somewhat  more  general  expression  for  probability  of  damage  as 


a function  of  weapon  usage « Is  given  by  Qallano  and  Everett. ‘ 

A comparison  of  the  results  of  this  formula  with  the  results 
of  calculations  against  Individual  target  elements  Is  given  by 
Schmidt*  for  a variety  of  oases.  A surprising  degree  of  agree- 
ment of  the  two  methods  Is  evident. 

Despite  the  venerablllty  of  the  Square  Root  Damage  Law,  no 
efforts  have  been  made  to  apply  this  general  approach  to  other 
situations.  That  Is  precisely  the  Intent  of  the  present  paper. 
In  this  attempt,  the  algebraic  structure  Is  presented  In  some 
detail  to  assist  In  understanding  the  nature  of  the  results. 

Chapter  II  presents  a derivation  of  the  damage  laws  when 
the  population  Is  In  a circular  normal  distribution,  or  uni- 
formly distributed  In  a disk  or  in  a ring.  Results  are  obtained 
from  situations  where  the  targeting  Is  optimal  against  one  type 
of  target  value  distribution  or  hardness,  but  the  actual  hard- 
ness is  of  a different  type.  Such  results  are  of  Interest  when 
two  separate  types  of  target  values  are  collocated,  e.g.,  where 
the  targeting  Is  optimized  against  the  Industry  In  a target 
complex  and  the  evaluation  Is  against  the  population. 

Following  this,  a more  generalized  derivation  of  the  damage 
law  allows  damage  calculations  for  value  distributions  that  are 
a function  only  of  radial  distance,  but  are  arbitrary  in  shape. 
Several  applications  of  th5 s derivation  are  presented.  Finally, 
some  factors  influencing  the  shape  of  the  damage  functions  are 
presented. 


^Robert  F.  Gallano  and  Hugh  Everett,  III,  Pe/iense  Modale  JV,  Family  of  Damaga 
FunoHona  for  Hultipla  Vaapon  Attaoka^  Lamda  Corp.,  Paper  6,  March  1967. 

*Leo  A.  Schnidt,  A Sanaitiaity  Analyaia  of  Urban  Blaat  Fatality  Calaulattana^ 
Institute  for  Defense  Analyses,  IDA  P-762,  January  1971. 
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Chapter  II 
ANALYSES 


A.  SQUARE  ROOT  DAMAGE  LAW  DERIVATION 

This  section  will  rederlve  the  Square  Root  Damage  Law  In 
some  detail  as  a basic  algebraic  theme  from  which  variations 
will  be  taken.  In  so  doing,  the  algebraic  structure  rather 
than  the  logical  basis  Is  of  prime  Interest. 

Assume  a population,  V^j,  Is  a circular  normal  distribution 
so  that  the  population  density  at  any  distance  r from  the  center 
of  the  distribution  Is 

V(r)  « VQ/2Tra^  exp(-r^/2a^)  . (1) 

In  subsequent  discus£"ions,  the  distribution  of  the  targeted 
vsxue,  Vip,  will  usually  be  different  than  the  actual  one;  in 
describing  the  former,  a will  be  replaced  by  in  the  latter, 
a will  be  replaced  by  o. 

Assume  a v;eapon  density  w causes  a fraction  of  fatalities 
P(w)  and  that  the  functional  form  ls‘ 


F(ta)) 

- 1 - e-‘'“  . 

(2) 

The  damage  at  any  point  is 

h(r,w)  ■ 

3 

U 

> 

(3) 

Through  the  introduction  of  a Lagrange  Multiplier  X,  one  can 


‘RiT  a discussion  of  the  relation  of  a weapon  density  to  a discrete  set  of 
targeted  weapons,  and  how  to  determine  values  of  the  constant  k,  see 
Schmidt,  op.  oit. 
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The  weapon  usage  W is  now  found  from 


r ® 

W - y l/k  ln(l/Y  exp(-r^/2a^) ) 27Tr  dr  , 


where  r„  solves 
c 


exp(-rQ/2a^)  - y • 


(12) 


This  is  readily  found  to  be 


W - 2iT/k  [ln(l/y)  • r?/2  - rJ/Ba^]  , 

o c 


W - 27r/k  Clnd/y)  • rV4]  , 

c 

W - (ir/k)a^  In^d/y)  , 


(13) 


using  the  definition  of  r twice.  A normalized  weapon  usage  is 
Introduced  by 


from  which 


X = kW/ira'"  , 


X » In'^d/y) 


y = exp(-  /X  ) . 


(IM) 


(15) 


A critical  point  in  achieving  a final  expression  for  survivors 
directly  in  X,  rather  than  two  expressions  in  the  parameter  y, 
is  solving  for  y in  terms  of  X. 

The  total  fatalities  are  found  from 


H * y* ' (1  - y exp(r^/2a^))  exp(-r^/2a^)  2iTr  dr  . 


Tnis  Is  1‘eadlly  Integrated  to 

H - Vq  (1  - exp(-r2/2a2))  - yVq  rl/2a^  . 

The  fatalities  are  expressed  as  a fraction  and  the  definition 
of  r.  Is  used  to  give 

H/Vq  - 1 - y - Y ln(l/Y) 

H/Vq  - 1 - Yd  + ln(l/Y))  . (16) 

Now  using  the  expression  for  y found  from  Integrating 
the  weapon  usage, 

HAq  - 1 - exp(-  ^ )(1  + /r*)  . (17) 

If  S Is  the  fraction  of  survivors,  then 

S - exp(-  /ST  )(1  + /ST  ) . (IG) 

B.  SURVIVORS  FROM  POPULATION  IN  A UNIFORM  DISK  TAR6ETE0 
AS  A GAUSSIAN  POPULATION 


TARGETED 


I 


Assume  a population  Vq,  for  which  damage  is  to  be 
assessed.  Is  uniformly  distributed  In  a disk  of  radius  R with 
population  density 
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r < R 


V(r) 


Vq/ttR^. 

0. 


* 


r > R . 


(19) 


The  targeting  Is  assuiaed  to  be  against  a Qausslan  dlstrl* 
buted  poDulatlon  given  by 

V - V^/2irB^  exp(*.r^/2B2)  (20) 

with  a damage  function  as  defined  In  the  square  root  law  deriva- 
tion. 

To  relate  the  assessed  population  to  the  targeting,  call 

* • R^/2B^  . (21) 

Since  the  targeting  Is  Identical  to  the  previous  section,  then 

X » In^dA), 

using  the  derivation  of  the  previous  section  with  a replaced  by 


To  compute  damage,  two  cases  mu&t  be  considered,  namely 

t 

Case  I;  r„  < R 

Case  II:  r^  > R 

c 

where  r , as  before.  Is  the  limiting  radius  for  u > 0.  For  Case 
I,  using  Equation  (11) 

h 2 2 2 

H ^ / Vq/ttR'^  • (1  - Y exp(r‘^/2$‘^))  2nr  dr  , 

0 


thus 


K/Vq  - 2/F^  [r^/2  - 0^'Yi-s-{,t>(r^/2B^)  - 1)]  . 


r!r.?!'^..,j..i,i|i..jli^l.i;PJ^^|i»^^,^,^ii..  J 


¥ 


|3. 

(A  ■ 

I; 


b 

1: 


m 


I 


ri 

c V 

■» 


Mow  using  the  definition  of 

H/Vq  • 26^/R^  Cln(l/Y)“  y/y  ♦ y}  - 1/^  Cln(l/Y)+  y - 1]  . 
In  terins  of  X,  using  Equation  l'15} 

H/Vq  - 1/#  [exp(-  /T  ) + /T  - 1]  . 

For  Case  II 

R 

H « y*  Vq/ttR^  * (1  - Y exp(r^/2B^))  2ffr  dr  . 


Thus 


So 


or 


H/Vn  - 2/R^  [R^/2  -0^  y exp(R^/2a2)  + B^y]  • 


H/Vq  “ !/<(•  [^  - Y exp(^)  + y3  > 


(22) 


(23) 


iH/Vg  - l/ifi  [.y>  + exp(-  /X  )(1  - exp(w))]  . (24) 

i 

To  dlstiriigulsh  between  Case  I and  Case  II,  we  notice  that 

! 

r.  • R at  the  Iborder  between  the  two  cases;  then,  by  definition 

C 

of 


exp(-n^/2B^)  • Y 

from  which,  using  Equation  (15)»  at  a critical  X 

? 

h ' 

Summarizing,  the'  fraction  fatalities  are  given  by 


(25) 


H/V, 


l/\p  [expC  - / X ) + /~3r*-  1], 


X < 


l/i»  i>fi  + exp(-  .nr  )(1  - exp(^))],  X . (26) 


i;| 

M f 


*, Ea^SW' " -S 
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In  tei^ms  of  survivors  t 

1 + IM  - l/^(exp(-  + /TT)  t X 

exp(-  »'T”)(cxp(#)  - X . (27) 

The  fraction  of  survivors  as  a function  of  X Is  presented 
In  Figure  1 for  various  values  of^.  The  dashed  line  separates 
the  two  regions  where  the  attack  radius  Is  less  than  or  greater 
than  the  target  radius  of  tht^  assessed  population.  Large  val- 
ues of w have  more  survivors  at  constant  X,  In  part  simply 
because  the  attack  Is  the  same  and  the  target  Is  larger. 

To  express  survivors  In  terms  of  weapon  usage  normalized 
to  the  assessed  rather  than  the  targeted  population,  let o be 
the  standard  deviation  of  the  targeted  population  about  the 
targeted  population  disk's  center;  i.e.. 


a 


2 


2Trr  dr 
irR^ 


SO 

o » r//T  . (28) 

2 2 

If  w Is  expressed  in  terms  of  o,  we  have  w ■ a /8  . To  normalize 
in  terms  of  the  actual  target,  let 

y » kW/iTO^  , (29 ) 

then 

Y - X/w  . 

Survivors  can  be  expressed  by  replacing  X in  the  previous 
formula  by  wY,  i.e., 

1 + 1/w  - 1/w  (exp(-  .^  + /^  ),  Y < w 

S = 

exp(-  /?Y  )(exD(W-  l))/y»,  Y > W.(30) 
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NORMALIZED  BY  TAR6ETIN6— ^ * R^/2B^.  X « kH/ir0 


Survivors  as  a function  of  Y are  presented  In  Figure  2. 
Since  the  normalization  Is  by  assessed  population*  the  lower 
curves  occur  with  more  efficient  targeting.  It  Is  Interesting 
to  notice  that  for  * 6*  where  the  assessed  population  Is  much 
more  dispersed  than  target  population*  the  targeting  Is  more 
efficient  for  heavy  attacks  (Y  ■ 10)  than  when  1 with  the 
two  populations  matched.  This  occurs  for  heavy  attacks*  since 
many  weapons  are  expended  outside  the  assessed  .copulation  disk 
where  ■ 1.  An  optimal  attack  against  a disk  has  a uniform 
weapon  density.  The  survivors  are  readily  seen  to  be  (see(35))5 

S = e“^/^  . (31) 

The  curve  labeled  'optimal  targeting*  in  Figure  2 exhibits 
this  formula  and  is  a lower  bound  for  any  of  the  other  curves. 

C.  SURVIVORS  FROM  A GAUSSIAN  POPULATION  TARGETED  AS  A 
UNIFORM  DISK 


Assume  a population*  Vq*  Is  Qausslan  distributed  with 
standard  deviation  o 

V(r)  - Vq/2tto^  exp(-r^/2o^)  . (32) 
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Figure  2.  FRACTION  SURVIVORS  FOR  POPULATION  IN  DISK  FOR  GAUSSIAN 

TARGETING  FOR  VARIOUS  VALUES  OF  WITH  WEAPONS  NORNALIZEO 
BY  POPULATION--^^  * r2/2b2  » Y » kU/na^ 


The  targeting  is  assumed  to  be  against  a disk  of  radius  R with 

IV^nR^  r < R , 

0 r > R . (33) 

Let  B be  the  standard  deviation  of  the  disk  about  the  center  as 
In  the  last  section,  so 

e • R/  /T  . 

p 2 

As  before,  define  ^ » a /B  to  describe  the  relative  size  of 
the  two  populations.  It  Is  clear  that  the  optimum  targeting 
against  a uniform  value.  Is  a constant  weapon  density,  thus 

W ■ ttR^u)  , 


then 

F(w)  « 1 - exp(-kW/-iTR^)  . 


Define  X by 


X = kW/irB^  , 


(34) 


then 


F = 1 - exp(-X/2)  . (35) 

The  damace  H Is  given  by 

R 'i 

H (1  - exp(“X/2))  Vq/2ito^  exp(-r^/2a^)  2irr  dr  . 

0 


This  Is  readily  seen  to  give 

H/Vq  • (1  - exp(-X/2))  (1  - exp(-P^/2o^))  , 


or 


H/Vq  - [1  - exp(-X/2)]  [1  - exp(-l/w)]  . (36) 
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For  survivors 


S ■ exp(-X/2)  Cl  - exp(-lM]  sxpC-l/^^)  . (37) 

To  express  survivors  In  ter'ms  of  the  assessed  population t 
we  define 

Y - kWAo^  . (38) 

Then 


Y - XA 


(39) 


so 

S ■ exp(“YW/2)  [1  - exp(>-lA)]  + exp(-lA)  . (^10) 

The  fraction  of  survivors  in  Figure  3 is  for  weapons  nor- 
malized by  targeting  and  in  Figure  ^ for  weapons  normalized  by 
population;  also  Figure  4 shows  the  Square  Root  Damage  Law 
which  represents  the  optimal  attack.  As  such>  it  must  be  a 
lower  bound  for  other  attacks.  As  can  be  seen  for  w ■ i/2,  the 
actual  survivors  are  close  to  this  lower  bound  for  almost  all  Y. 


D.  SURVIVORS  FROM  A GAUSSIAN  POPULATION  TARGETED  AS  A 
GAUSSIAN  POPULATION  OF  A DIFFERENT  DISPERSION 


Assume  a population,  Vq,  Is  Gaussian  distributed  with 
standard  deviation  a so 
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iir'*ni  i ■jSfSir  * '"mi'ffi 


" ' . »'  i,j.i,.iiiini....ing- 


Figure  3.  FRACTION  SURVIVORS  FOR  GAUSSIAN  POPULATION  FOR  TARGETING  OF 
UNIFORM  DISKS  OF  RADIUS  R WITH  WEAPONS  NORMALIZED  BY  TARGET 
= 2ir2/R2  * «r2/02^  X = kW/irB^ 


i/y3  . 


(45) 


and  recalling  the  definition  of  r.  gives 

w 

H/Vq  - (1  - “ Y/(1  -^)[1  - 

2 

Now  defining  X as  kW/TrP  , and  substituting  for  y using 
Equation  (15), 

H/Vq  = 1 - exp(-  /T/*) 

-[exp(“  /r~)/{l  -^)]Cl  - exp(-  /r'/^)/exp(-  /JT")]  . 

(46) 


For  survivors, 

S = (^/^-l)  exp(-  vOT'/v)  - exp(-  /3r'/('^-l))  . (47) 

To  normalise  by  the  actual  population,  let 

Y = kWAa^  . (48) 


So 

Y = XA  . (49) 

Then 

S = ('P/w-l)  exp(-/YA  ) - l/(«^-l)  exp(-  ✓ Y<^  ) . (50) 

It  Is  Interesting  to  notice  that  If ^ Is  replaced  by  1 A , 
the  same  formula  is  obtained.  This  Implies  that  if  the  target 
population  is  either  too  small  by  a factor  of  X or  too  large 
by  a factor  of  X,  relative  to  the  actual  population,  the  same 
number  of  survivors  Is  obtained.  The  fraction  of  survivors  is 
presented  In  Figure  5 with  weapons  normalized  for  targeting  and 
in  Figure  6 with  weapons  normalized  for  assessed  population. 

E.  SQUARE  ROOT  INJURY  LAW 

Assume  a population  is  Gaussian  distributed  so  that 

V(r)  = VQ/2-na^  exp(-r^/2o^)  . (51) 

Assume  the  population  is  targeted  so  that 

F(a))  = 1 - . (52) 
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Igure  5.  FRACTION  SURVIVORS  OF  GAUSSIAN  POPULATION  WITH  STANDARD  DEVIATIOi 
WITH  TARGETING  AGAINST  GAUSSIAN  POPULATION  WRH  STANDARD  DEVIATE 
WITH  WEAPONS  NORMALIZED  FOR  POPULATION  0^/02,  % = kH/ir02  / 


ppr.  _ 

•'  ¥■ 

■ ^'.i 

t 

V 


j^rrg'TOjwiiaiiiaffmH 


Now  consider  another  eya’^'uat Ion  function j 


F^Co))  = 1 - e 


-k^w. 


This  function  will  he  us 3fi 'for  two  casesj' 

Case  I:  The  actuuj  vulnerability  is 'k^  and  the  effect 

of  the  attacker  using  an  3<hcorrect  vulnerability 
' ■ i^^  ,pf  lnt\‘rest  j ,, 

Case  ,11:  IJThe  k^ '''^,^'>r.esents  a casualtj^' rate  and  we  wish  > 
to  deter, iffline  injuries  (casualties  minus  fatal- 
ities). 


In  either  case,  define 


. 


C = k^/k 


(54) 


For  fatalities,'  when!  P(u>)  is  used  directly  with  k,  the 
square  root  law  derivation  is  used,  leading  to 

H/Vq  = 1 - exp(-  t^T  )(1  + /T  ) , 


where  X Is  defined  as  kW/ira  . 

Call  H the  \3i}.ue  obtained  using  k . 


H. 


r ° 


(55) 


Figure  7.  FRACTION  SURVIVORS  AS  A FRACTION  OF  WEAPON  USAGE  FOR  VULNERABILITY 
INDEX  AND  TARGET  k WITH  WEAPONS  NORMALIZED  FOR  ATTACK— C = k./k 


Then 


Sq  - ^ 5 i «*p("  «*p(-  ) • (63) 

Replacing  £ by  1/C  In  this  formula  leaves  survivors  unaf- 
fected} which  Implies  that  either  an  underestimation  or  an 
overestlmatlon  of  vulnerability  will  give  the  same  result. 
Moreover,  this  formula  Is  Identical  to  that  for  survivors  of 
two  Gaussian  populations  of  different  dispersions  In  the  pre- 
vious section' with  ^ replaced  by  C*  The  curve  of  survivors  In 
that  case.  Figure  6,  also  applies  here. 

For  Case, II,  call  Injuries  the  difference  between  fatal- 
ities and  casualties,  so 

I - H^/V„  - H/Vg<  ^ (6H) 

Then  since  H/Vq  **  1 - y - Y ln(l/Y),  we  h^e,  from  Equation  (59), 

I = + Y ln(l/Y).  (65) 

Thus  I 

I - -.■ex£.(-A>nC  ) + ^ exp(-/T  ) . 

or 

I - exp(-  /T.)  ) + /x  . (66) 

Injuries  are  presented  In  Figure  8 for  values  of  C greater 
than  1.  In  this  Figure,  the  attacker  Is  attempting  to  maximize 
fatalities.  If  t^e  attacker  is  attempting  to  maximize  total 
casualties — Injuries  plus  fatalities — then  valuesfof  C are  less 

r 

than  Ij  l.e.,  the  vulnerability  of  the  evaluated  population  is 
less  than  assumed  by  the  attacker.  The  Injuries  In  this  case 
are  presented  In  Figure  9. 
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Figure  9.  INJURIES  AS  A FUNCTION  OF  WEAPONS  OSAGE  WHEN  ATTACKING  FOR 
TOTAL  CASUALTIES--^  - k^/k.  X « kV/na^ 


SURVIVORS  FROM  A GAUSSIAN  POPULATION  TARGETED  AS  A UNIFORM 
POPULATION  IN  AN  ANNULUS 


ASSESSED 


Ml-Tt-U 


TARGETED 


Assume  a population,  Vq,  is  Gaussian  distributed  with 
standard  deviation  o so 

V(r)  « Vq/2ito^  exp(-rV2a^)  . 

The  targeting  is  assumed  to  be  against  a uniform  ring  with 
inner  radius  L and  outer  radius  H.  In  order  to  normalize  the 
evaluation  in  a way  comparable  to  other  situations,  a “center 
of  gravity”  and  "standard  deviation"  of  the  ring  is  defined 
by 


H 

m » 1/A  J'  r • 2TTr  dr  , 


(67) 


2 2 

where  A is  the  ring  area  = iv(H  - L ),  and 


» 1 A y*  (r  - m)^  2TTr  dr  . 


(68) 


- • ir' 


Thli  yields 


_ _ 2 - L* 

in  • X N » 

3 - L‘ 


„2  _ 1 - L 

® “ 5 ::5 — n 


(69) 


(70) 


- H'-  - L‘ 

The  expression  for  B Is  a measure  of  dispersion  about  the 
"center  strip"  of  an  annulus.  Thus^  as  the  Inner  radius  of 
the  annulus  approaches  0,  B does  not  approach  the  value  com- 
puted for  a disk  which  Is  about  the  center.  In  fact,  with 
L ■ 0,  we  have 


®ring  “ ^*^3  ®(aisic 


There  Is  no  simple  solution  for  H and  L in  terms  of  m and 
B,  but  a simple  Newton's  method  type  of  procedure  will  allow 
for  quickly  determining  numerical  values.  Figure  10  shows  m 
and  B as  a function  of  H and  L.  As  the  thickness  of  the  ring 
becomes  small  compared  to  a radius,  the  ring  curvature  Is 
insignificant  in  determining  m and  B,  so  that  we  have  approxi- 
mately 


H « m + , 

L w m - , 


B/m  « 1 , 
B/m  <<  1 . 


The  departures  of  the  contours  of  constant  m and  B in 
Figure  10  from  straight  lines  where  L/H  is  small  indicate 
departures  from  this  approximation. 

The  ratio  of  target  to  actual  standard  deviation,  B/o, 

we  shall  call  p,  and  the  ratio  of  ring  center  to  standard 

deviation  m/o,  we  shall  call  n*  At  small  L,  the  w used  previ- 

2 

ously  is  related  to  p by  w ■ l/9p  if  6 is  substituted  for 
3B. 


The  total  weapon  usage  W is  given  by  weapon  density  times 
ring  area,  or 


.1  aWfcui  zaial  kua 
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W - - L^)to»  . 


Call  X the  weapon  usage  normalised  to  the  target,  or 

X « kWAB^  • k«(H^  - . (71) 

Call  Y the  weapon  usage  normalised  to  the  assessed  population, 
or 

y ■ kWAB^  • kta>(H^  - l?)/o^  . (72) 

The  succeeding  calculations  will  be  In  terms  ox''  Y but  could 
be  readily  converted  to  X by  changing  o to  B.  The  calculations 
will  be  done  explicitly  In  terms  of  H and  L,  which  are  functions 
of  m and  B.  As  the  expression  for  Y shows,  these  functions  can 
be  considered  as  normalised  by  c;  l.e.,  the  expressions  are  In 
terms  of  HA,  LA,  P and  n* 

The  fraction  survivors  k Is  given  by 


P ■ 1 - exp(-kw)  ■ 1 - exp(- 

The  value  destroyed  is  given  by 
H 

H - / F . V dA  i 


- L^A^ 


•)  . 


(73) 


H 


• f (1  - expC ^ ^ — 2 — g-])  Vq/2ttct^  exp(-r^/2a^)2Trr  dr 

*>  H ^ - L /cr 


Integrating,  and  letting  the  fraction  of  survivors,  S 

5 2 


1 - H/Vq,  gives 


S • 1 - [1  “ exp 


(- 


(HA  ) - (LA) 


[exp  (-  “ exp  (-  -^)]. 

^ / 2o^  2d*^ 


(7^) 


(A  similar  formula  based  on  S is  attained  by  substituting  X for 
Y and  B for  a. ) 
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Fraction  survivors  are  Indicated  In  Figure  11  for  the  annu- 
lus center,  m,  at  a constant  value  with  varying  annulus  width, 
Bi  and  In  Figure  12  for  the  annulus  width,  B,  cor.  ttant  and  vary 
Ing  distance  of  the  center  m.  In  these  figures  a curve  for 
L > 0 Is  given  which  reduces  to  the  disk  calculations  presented 
earlier. 

When  L/H  Is  near  1 and  the  approximate  formulas  for  H and 
L can  be  used,  we  have 

S « 1 - [1  - exp(-Y/(l»  /Inp)  )]Cexp(-(n  - /5p)V2) 


- exp(-(ri  + /5p)  /2)] 


(75) 


G.  SURVIVORS  FROM  A POPULATION  IN  AN  ANNULUS  TARGETED  AS  A 
GAUSSIAN  POPULATION 


fARGETED 


ASSESSED 


i-n-Tt-iT 


Assume  a population  is  of  constant  density  In  a ring  of 
inner  diameter  L and  outer  diameter  H.  Let  m and  o be  the 
mean  ring  distance  and  o the  standard  deviation.  The  relations 
between  H and  L and  m and  o are  the  same  as  In  the  previous 
section,  with  o replaced  by  B.  The  targeting  Is  assumed  against 
a Gaussian  target  distributed 

V » exp(-r^/2S^)  . 

As  In  the  square  root  law  derivation,  the  weapon  usage  Is 


WITH  CONSTANT  WIDTH— p « B/a  * 0.25 


with 


X - kW/irP^  - ln^(l/Y)  , 

Y - air&^x/kV-Q  . 

The  fraction  damage  P is 

P “ 1 - Y exp(r^/26^)  , 

and  extends  to  a distance  where  solves 

c c 

Y * exp(-r|/2B^)  . 

Por  r^  less  than  L there  is  no  damage  and  the  fraction 
survivors  is  1.  Otherwise*  for  damage 

J 

H - y [1  - Y exp(r^/2S^)]  CVqACH^  - L^)]  2xr  dr  , (76) 

L 

Where 

J - min(r,,H)  . (77) 

This  readily  gives 

S - 1 - H/Vq  - 1 - - 2 exp  [exp(J^/2B^) 


- exp(L^/20^)]}  . 

If  r^  > H this  simplifies  to 


3 » 2 exp  texp(H^/2B^)  - exp(L^/2B^)  J . 

- L 


(78) 


Por  r„  < H this  can  be  written 
c 

S - 1 - g - 2 exp(-  /T  ) • 0^  • exp(L^/20^). 

H - L H - L 


(79) 


3^ 


— a.Uia.Li 
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In  terms  of  the  assessed  population,  the  normalization  becomes 

Y - kW/irB^  - x/i()  (80) 

where  « B^/B^.  The  fraction  of  survivors  as  a function  of  Y 
is  shown  In  Figure  13  for  an  anrtulus  of  constant  ratio  of  outer 
to  inner  diameter  of  2.9.  When  the  targeted  population  is  not 
dispersed  enough,  the  weapons  are  concentrated  in  the  hole  in 
the  annulus.  As  the  standard  deviation  of  the  targeted  popu- 
lation is  Increased,  the  tarpeving  becv>ines  closer  to  optimal 
targeting  until  finally  the  v,.,apons  are  too  dispersed  and  thi 
targeting  effectiveness  drops  off. 

H.  AN  ALTERNATIVE  SQUARE  ROOT  LAW  DERIVATION 

In  this  derivation,  the  geometric  context  ,of  the  square 
root  law  will  be  emphasized.  As  before,  let 

V(r)  , - V^/2Tra^  exp(-r^/2o^)  . 


Let 


so 


Then 


» 2a^ 


/2^  . 


(81) 


V(r)  « exp(-r^/6^) 


-1 


i 1 
II 


I 


Now  call 


.. 

. :i 


1 . ^ 


Tffi' 


(82) 


so  A,j,  is  the  circle  which  contains  1 - 1/e,  or  63  percent  of  the 
population. 


Let 


A “ wr 
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SURVIVORS  FROM  A UNIFORM  ANNULUS  POPULATION  TARGETED  AS 
GAUSSIAN  POPULATION  FOR  H/L  =2.9 


Then 


V(r)  - Vq/A^  exp(-A/A^)  . (83) 

As  before,  let 

P(w)  - 1 - e“*^“  . 

The  optimal  allocation  yields  the  requirement 

- X/kV,  for  X/kV  < 1 . 

Now  let  s,(u)  denote  the  fraction  of  survivors,  locally,  so 

s(w)  » e"^'^  - X/kV  . 

Call  G « sV  the  survivors  per  unit  area.  Then 

G = V • X/kV  = X/k  . (84) 

Thus  the  survivors  per  unit  area  are  a constant  for  the  area 
attacked. 

The  attack  continues  to  a radius  where  X/kV  * 1.  Call  the 
area  enclosed  A^.  At  this  radius  the  target  value  equals  the 
survivors  per  unit  area.  Thus, 

V/G  = 1 , 

so 

Vq/GAij  exp(-Af,/A,p)  = 1, 


or 

ln(VQ/GA^)  = Ap/A^  , 

Define,  as  before,  a dimensionless  weapon  usage  by 

^C 

X = kW/ira^  = 2/A,p  J"  kw  dA  . 

0 


Now,  since 

» G/V  , 


(85) 


(86) 
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Thus 


ku>  - ln(V/a)  - ln(VQ/aA^)  - A/A^  . 


c 

X - 2/A^  f (ln(VQ/QA,p)  - A/A^)  dA 


2/A^  ln(VQ/aA^)A^,  - 2/A^^  A^^/2 


2 Ag/A^  (ln(VQ/GA^)  - 1/2  A^/A^)  . 


Since  A^/A^  - MV^/QA^)  . 


or  alternatively. 


X - (A^/A^)*^  . 


X - ln'^(VQ/GA^)  . 


(87) 


Divide  survivors  Into  the  area  In  the  circle  A^,  and  the 


rest  of  the  plane 


S ■ + S2  • 


Si  ■ GA^  . 


CO  OD 

^2  “ f 1 • V dA  - y Vq/A,p  exp(-A/A^)  dA  . 


Vq  exp(-A^/A^) . 

S » GAp  + Vq  exp(-AQ/A,p) . 


Normalizing, 


S/Vq  - GAq/Vq  + expC-Ap/Arp)  . 

Now  S will  be  computed  both  in  terms  of  G and  Aq. 


(88) 
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First,  geometrically: 


Now, 


S/Vq  ■ QArp/VQ  • Aq/A^  + exp(-Aj,/Aj)  . 

In(QA^VQ)  • “Ag/A^  . 


QA^/Vq  ■ expC^A^/A^)  , 
S/Vq  - expC-Ag/Aij,)  Aq/A^  ^ exp(-A(,/A^)  , 
" ® XP  ( “Aq/A^  ) ( Aq/A|J| 

And  since 

Ag/A^  « «/7“  , 

S/Vq  » exp(-/5^  )(»^  +1)  . 

Now  In  terms  of  survivors 


- GA^/Vq  ^ln(VQ/GA^)  + l)  . 


Now 


GAij/Vq  » exp  (-/ST  ) , 

so  * expC-fOT”)  (/T”  +1)  . 


(89) 


(90) 


1.  A GENE[ A.iZATION  OF  THE  SQUARE  ROOT  DAMAGE  LAW 

Suppose  value,  as  a function  of  area  covered,  is  an  arbi- 
trary decreasing  function.  Then  the  use  of  weapon  density  can 
still  be  used  for  a more  general  form  of  the  square  root  law. 


ao 

6-*'“  . for  5^  < 1 . 

Call  a the  local  fraction  of  aurvivora  and  Q the  local  number 
of  aurvivora  per  unit  area.  Then 

G - Va  - V(l-f)  - Ve‘^“  - V-^  » ^ . 

Thud  Q la  a constant  In  the  area  attacked.  The  attack  la  made 
out  to  a radius  where 


or 

Call  the  area  attacked 


The  weapon  usage  W la  given  by 


kV  • ^ 


Then 


(5r) 


3A 


(96) 


w 


bidA 


Let  X be  the  normalized  weapon  usage,  as  before,  so 


ill 


(97) 


■' 

^ ^ ln(f(x))dx  . 

0 

(98 

Now  define 

h(x)  - J 

r* 

' ln(f(t))dt  . 

0 

(99: 

Then 

X • 2^1n 
At 

fe)  * '"ft)  • 

(10( 

The  total  survivors 

Is  the  sum 

1 S,  of  the  survivors  in  the 

area 

attacked  and  In 

the  unattacked  area  with 

Tiiie  firttctlon  of  survivors  is  then 


S 


55 


+ 1 


J^o\ 

" *lzri  • 

V*T  / 


(101) 


J. 


AN  APPLICATION  OF  THE  GENERALIZED  DAMAGE  LAW  TO  VALUE 
A LINEAR  FUNCTION  OF  AREA 


Let  value  be  a linear  function  of  area,  say  of  the  form 

^ ^ - 4)  • * - 

0 » A ^ Ap  • 


V(A) 


, 2 
'0  r;i 


(102) 


For  a circularly  symmetric  value  function,  value  would 
decrease  quadratically  with  distance.  After  normalizing  with 
the  area  A^,  we  have 

V, 


4 *4) 


iL.  < i 

Aj  - § 


0 , ' ifK/ 

s 

where  e ■ 1 -/l/e  a 0, 393^6'93i|. 


‘T 


> Y 


(103) 


(lOH) 


Then 


g(x) 


-2  2 

2ex  -ex 


X < i 


* ' f > 


and  it  is  readily  seen  that  g(l)  »»  1 - as  desired  for 

w 

normalization.  Moreover,  in  terms  of  t,he  original  triangle. 


A^  — eAp  . 


H3 


For  h(x)  we  have 


h(x)  ■ ^x  - log (28  - 28^x)  “ + g log(2S)  . (1Q5) 

Now  suppose  a value  of  Q (survlvors/unlt  area)  is  given. 
Then  with  defined  as  the  area  under  attack,  call  x^  ■ A^/A,p. 


Then 


1 

e 


QA, 


T _1 
0 28 


Now  from 


2x^  In 
c 


(qa°)*  » 


we  have 


(106) 


and  from 


IT  " ^ 
''o 


+ 1 - g(x.)  , 


we  have 


''o 


1 ®^T  1 /^^t\ 


(107) 


Numerically,  this  gives 
X » 5.0829881  In 


s/Vq  » 2.5^149^0 


6.i.5919n(^)-  6. 

/OA„\  /QA„\2 

1.6i4798o(-^)  . 


3008972 


In  this  case  no  explicit  function  ^(x)  is  presented. 

0 , GA„ 

Solutions  must  be  obtained  in  terms  of  the  parameter 
The  valid  range  of  area  attacked  is 

0 1 Xq  - I * 


nn 


1 

vr 

e 


which  gives 


QAm 

2g  * 0.79693868  > ^ > 0 . 

- Vq  - 


The  survivors,  as  a function  of  attack  Intensity,  Is  pre- 
sented in  Figure  The  square  root  law,  with  Qausslan  tar- 

gets, Is  Included  for  comparison.  As  can  be  seen,  practically 
no  difference  occurs  between  the  two  curves  until  the  survi- 
vors are  under  50  percent.  Then  the  larger  amounts  of  area  at 
low  value  density  in  the  Gaussian  targets  begin  to  signifi- 
cantly affect  the  attack  effectiveness.  In  Figure  15,  the 
survivors,  attack  size,  and  several  other  parameters  are  pre- 
sented as  a function  of  the  area  attacked.  As  is  clear  in  the 
Figure,  not  until  almost  all  the  area  is  attacked  does  the 
weapons  usage  become  relatively  large. 


K,  AN  APPLICATION  OF  THE  GENERALIZED  DAMAGE  LAW  TO 
VALUE  AN  EXPONENTIAL  FUNCTION  OF  AREA  TO  A POWER 

Let  the  value  as  a function  of  area  be  given  by 

V(A)  - VqH  exp{-(6A)’^}  , 0 < n < « . (108) 

2 

The  area  A is  Ttr  , so  this  can  be  written 


I 


V(A)  = exp^-{BTTr^}  ^ . 

For  n = 1 the  standard  Gaussian  shaped  curve  is  obtained. 
For  n « ^ an  exponential  decay  is  obtained,  a and  B are 
parameters  whose  values  will  be  determined  by  normalizing  the 
value  function.  The  damage  will  be  calculated  using  the  second 
altei*native  derivation.  To  normalize,  call 


Then 


, f(x)  = a exp(-(Bx)”)  , ' 

w.iere  now  a and  ,6  ,ca,i:i  be  expressed  as  functions  of  n.  Now  for 
valu^' we  can  also  write  ' 


V(A)  - I . 


(117) 


and  preserve  normalization . Here  Am  is  the  area  within  which 
1 - — of  thr  population  resides.  Then.,x  » A/Am  becomes  area, 
expressed  in  uni^s  of  A^.  , . 

Now  let  a number  of  Survivors  per  unit  area,  G be  given.  We 
wish  to  find  ,an.  x^(»  A^^/A^i)  so  that 


GAm: 

f(x  ) » - X . 

, .•  '0  „ , , 


(118) 


Thus  x^  solves 
c 


a exp(-(‘0x  )'^)  , » I . 


(119; 


Since  > 0,  then  a condition  on  X*  is 
Q . 

X <_  a 

The  function  h(x)  is  defined  as' 


h(x)  « r ln(f(t))  dt  . 


(120) 


Here 


h(x) 


f 


ln(q  exp(.-(6t)”}  )dt 


X in  (a,)-  P-Vt 


(121) 


k9 


In  particular, 


Now 


■ V , V, 

r+T  • 


, X - ln(l/X)  ,+  2h(x^)  . 

K*  C I - ' , 

Finally,  for  survlvoi'^i  ^ ,: 

. ' • ■ , s X 1,+^  1 - :i(x) , . '■' 

; The  following;  table  compares  exact  , values  of 
values  fro»  the  fitted  curves.  ■ , 


(122) 


(123) 


fl24) 
8 and  with 


n 

*^tiiact 

®f1t 

— 

:1s)  ; , 

5n  ’ ? 

1.  -‘i' 

nrr 

“fit 

0.125 

27.100.000 

197,540 

5.040 

0.0014878 

-5  51044 

672.1230 

84.60 

0.25 

357.473 

201 

' ' • 5 

0.0<37|4 

-2.701 

14.89 

14.76 

0.5 

4.5224 

5.0 

1 

0.43258 

-0. 83777 

2.3112 

2.58 

0.75 

1.55201 

1.6! 

0.59290 

8.7572 

•2.5506 

1.3036 

1.18 

1.0 

1.00231 

0.995 

1.Q0 

0 9977 

•0.0C23 

1.0C23 

1.00 

1.5 

0.716534 

0.59;* 

1.i640 

1.2595 

0.23079 

.0.71939 

0.0529 

2.0  ' 

0.63775 

C.523 

1.7724 

1.38911 

0.3289 

0.7197 

0.7897 

4.0 

0.55720 

0.691 

3.4245 

1.5435 

0.4341 

0.6478 

0.7058 

a.o 

0,597153 

0.502 

7.534 

1.5771 

0.4556 

0 6341 

0.6706 

54.0 

0.5274T 

0.63^ 

63.44 

1.6790 

0.4574 

0.5329 

0.6437 

Figures  ^ 6 and  17  present  values  as  a function  of  distance 
for  n > 1 r.ud  n < Ij  that  is,  the  curves  present  f(x)  as  a func- 
tion of  /x.  The  Integral  of  the  value  g(x)  is  presented  in 
Figure  18  as  a function  of  distance  /x  for  the  extreme  values 
of  ri  considered.  In  Figure  19  the  fraction  survivors  as  a 


‘The  "exact”  vaiuea  of  8 and  a are  found  throtigii  solving  the  equation 

'■ " ' 1 /*6’ 

g(l)  I exp(-y‘')  dy  . 

® •'0 

This  was  acccnipli4,»hed  numerically  by  evaluating  the  Integral  by  Slnpson’s 
rule  iJisiiig  1,000  steps  for  differing  values  of  8 until  the  desired  value  was 
obtained.  Since  at  n“l,  8 aJ-jould  be  1 exactly,  the  quoted  value  gives  arv 
Ihdicetion  of  the  error  In  this  numerical  prcced’ure. 
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VALUE  AS  A FUNCTION  OF  DISTANCE  FOR  n GREATER  THAN 


VALUE  AS  A FUNCTION  OF  SQUARE  ROOT  OF  DISTANCE  FOR  n LESS  THAN  1 


function  of  normalized  weapon  usage  Is  presented.  For  n ■ 1, 
the  curve  Is  the  usual  square  root  damage  law.  For  n ■ 64, 
the  value  Is  almost  uniform  and  then  rapidly  drops  so  the  dam- 
age Is  close  to  that  of  a uniform  disk  optimally  targeted;  l.e., 
the  fraction  survivors  Is  an  exponential  function  of  weapon 
usage.  For  n ■ 1/2,  the  value  Is  an  exponential  function  of 
distance.  For  values  of  a less  than  1/4,  no  curves  are  pre- 
sented since  the  decay  of  value  with  distance  Is  extremely 
rapid  and  numerical  Integration  using  Simpson's  rule  over  Inter- 
vals (the  algorithm  Implemented  here  for  numerical  evaluations) 
becomes  tedious. 

L.  FACTORS  INFLUENCING  THE  PROBABILITY  OF  KILL  AS  A FUNCTION 

OF  WEAPON  DENSITY 

In  this  section  some  rational  for  the  use  of  a weapon  den- 
sity Is  presented.  This  Is  done  first  by  studying  the  survival 
probabilities  when  a larger  number  of  weapons  arrive  at  random 
in  a large  area.  Next  the  probability  of  survival  from  weapons 
delivered  In  regular  patterns  in  an  area  are  considered.  This 
Is  done  analytically  for  weapons  where  the  probability  of  kill 
for  the  weapon  Is  one  Inside  a weapon  radius  and  zero  outside, 
and  numerically  for  several  typical  types  of  assumed  weapon 
kill  probabilities  as  a function  of  distance. 

Suppose  weapons  are  delivered  In  an  Infinite  plane  so  that 
any  point  In  the  plane  is  as  likely  to  receive  a weapon  as  any 
other.  Then  in  any  specified  area  the  number  of  weapons  deliv- 
ered Is  a Poisson  process  with  the  distribution  of  number  of 
weapons  given  by 

P(„)  . 6-“*  . (125) 

where  A is  the  size  of  the  area.  Since  the  expected  number  of 
weapons  In  a unit  area  Is  to,  we  will  call  to  the  weapon  density. 
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Consider  any  monitor  point  In  the  plane  and  let  this  moni- 
tor point  be  surrounded  by  a set  of  N non- overlapping  but  touch- 

tiH 

Ing  annuli,  where  the  J annulus  is  centered  a distance  r^ 
from  the  monitor  point  and  has  a thickness  dr^ . 

The  area  of  the  annulus  la  given  by 
Aj  « iT(2r,drj  + dr^^)  . 

The  probability  of  n weapons  arriving  In  this  annulus  Is 


P(n) 


-tjAj  (to)A^  ) 


n: 


n 


The  overall  probability  of  kill  of  the  monitor  point  from 
weapons  In  the  annulus  Is  the  sum  of  probabilities  of  kill 
for  one  weapon  arriving,  two  weapons  arriving,  etc.,  or 


-wA, 


...  *(i  - 


(1  - P^(rj))» 


)(-,  ■ 


The  rings  form  a mutually  exclusive  set  of  areas  covering 
the  plane,  except  for  a distant  final  outer  section  surrounding 
the  rings.  Then  the  probability  of  survival  of  the  monitor 
point  is  given  by  the  product  of  the  probabilities  of  surviving 
the  weapon  arrival  in  each  ring. 

Now  let  N become  Indefinitely  large  so  (1)  for  all  J for 

J » 1,2,.,.N,  11m  A.  ■ 0 and  (2)  lim  r,j  « «.  It  Is  clear  that 

N-^ 

this  can  be  done  with  a countable  set  of  rings;  for  example, 

oyii 

let  Tj  - J/N-’  j « 1,2,...N.  Then  the  area  of  each  ring  is 


and  13m  A.  » 0. 
N-*-« 
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S • lljn  n U - ^ J 

J-1 

Nr  2 1 

- llm  n U - ■*■  ' 

N-®-  j.3^L 

Now  taking  logarithms  of  both  sides, 

Nr  2 1 

in  S - in  llm  n U - ^k^^j 
N-- 

- llm  £ ln(l  - O(drj^))  . 


2 3 

in  (1-x)  ■ - (x  •<•  ^ + ^ + . • • ) |xl  < 1 . 

1 

hus,  y 

N , 2\ 

In  S » llm  ^ -Pj^(rj  )w2TTrjdrj  + ) . 

’**“  3-1 

Now,  assuming  the  P?,(r^ ) are  such  that  the  sum  converges 
• ^ J 

which  lS'^;fcertalnly  true  for  all  practical  damage  functions. 


(126 


F! 


■ i 


11 


Now  J Pj((A)  dA  Is  often  equated  to  irWR^  where  WR  Is 

called  the  weapon  radius.  The  circular  area  Is  the  area  of  a 
cookie  cutter  damage  function  which  does  the  same  damage. 

Thus,  finally, 

3 e g-wir(WR)  ^2 

and  2 

P . 1 . . 

This  Is  the  form  of  the  damage  function  of, Equation  2 which 
underlies  the 'square  root  law  formulation. 

A suggestion  by  John  Donelson*  gives  the  same  result  In  a 

much  shorter  fashion  and  In  a way  which  makes  the  answer 

Intuitively  clear.  The  expected  fraction  of  the  plane  covered 

2 

by  lethal  effects  Is  wTr(WR)  . A small  change  In  this  fraction, 
l.e.,  AwttCWR)*^  gives  a proportional  fraction  change  In  the 
survivors . So 


or 

so 


- ^ - AwTr(WR)^  , 


-Ain  S ■ Au)ir(WR)*^  , 


S . ^-ot(WR)2  _ 


(128) 


The  random  pattern  of  Poisson  arrivals  just  discussed 
might  be  appropriate  for  an  attack  against  certain  types  of 
valuable  targets  in  an  area  where  population  happens  to  be  cal- 
culated, and  where  the  location  of  the  population  Is  not  sig- 
nificantly influenced  by  the  location  of  the  other  valuable 
targets.  On  the  other  hand.  If  an  attacker  were  Interested  In 

^Btonrfirly  with  Institute  for  Defense  Analyses,  currently  associated  vdth 
Science  Applications,  Incorporated. 
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i 
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attacking  an  entl:'e  urban  area,  then  he  might  place  his  weapons 
In  a regular  pattern.  Any  target  In  the  urban  area  Is  some- 
where within  the  pattern,  but  unless  some  reason  is  given  to 
fix  the  location  of  the  weapon  patcern,  the  target  may  be  any- 
where within  it . 


^ ! 


i i 


Suppose  weapons  are  placed  at  the  Intersections  of  a square 
grid  with  mesh  length  S..  Then  the  weapon  density  w is  simply 
(0  ■ 1/a  . If  the  probability  of  kill  Is  1 inside  some  weapon 
radius  WR,  and  0 outside,  then  the  expected  kill  is  simply 
determined  by  the  fraction  of  the  plane  covered  by  weapon  circles, 
Call  p ■ WI, the  fraction  coverage,  and  expected  kill  is 
readily  seen  to  be 


TfP 


- 2cos“^(i)^  , 


0 < p < 


1 < p < 


< 


1 

/2 


(129) 


The  fatalities  can  be  given  as  a function  of  the  fraction 
coverage  of  the  plane,  C,  to  be  consistent  with  the  previous 
section.  Then 


iT(VfR)^U) 


ir  2 
TfP 


(130) 


If  the  weapons  are  dropped  in  a more  closely  packed  pattern 
of  equilateral  triangles  of  side  length  I,  then  w » 2//3^  . 
Again,  calling  p « WR/S,/2  for  fraction  fatalities. 


'TK2 


0 < p < 1 


4(f  ■ , 


1 , 


The  coverage  is  given  by  C = 


1 < p < 


P 1 7Z 
/3 


/3 


IT  p2 

2/3 
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LINEAii  PLOT  OF  SURVIVORS  AS  A FUNCTION  OF  FRACTION  COVERAGE 
FOR  SEVERAL  TYPES  OF  ATTACK  ^ 


FRACTION  COVERAGE  OF  PLANE  (PERFECT  TARGETING) 


Figure  21.  SEMI-LORARITHMIC  PLOT  OF  SURVIVORS  AS  A FUNCTION 
OF  FRACTION  COVERAGE  FOR  SEVERAL  ATTACK  TYPES 
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' '"••'!  7,'  , / 


/r^'“fi'm,f;V‘' 


The  fraction  survivors  ei^^fujiction  ,pt  coverage  Is  pre- 
sented as,  a linear  ,plo't  1^  Mi^ure  20'a^d  as  ',a  isemi-legarlthmlc 
plot  in  figure  Also  presented  are  survivors  for,t|^e  raridOKi  ' 

delivery,  considered  eaplier,  T^d  bottom  straight  line  in 
Figure  20  represeiits  thfe  survivors; -if',  weapon  effects,  could  be  , 
perfectly  p'laoed>  e,g. , square  iefehal- effecta' ai^eas  covering  ' 
the  plane  for  weapdns  plape^,  in' '^hs  tenter' of  each  square.  As 
can  be  seen  * there  is  no  overlap  for'  the  equ'ilat eraT  triai^le, 
pattern -until  0.’9  -of  -the'.' plane,  -is  covei^di' ''b'y-''w'teapohs-..,effe;ots.^ 

Even  then,  there- is  less  than  four,  .percent  difference  between 
this  coverage  and  t,he  no  overlap  variation.  ' < 

If  the  fracticn  damage  as  a function  of  distance  is  a 
gradually  decreasing  function  of  the  distance,,  t^ien  effects  of 
overlap  must  be.  considered.  Figures  20iiand  21  frhiiijW  cailoulatlp,ns 

' ' ' * , . , ■ • ■ • I ■ j //  . ^ ,1 

for  four  eases;., of  Figure  22  with  weapon^, 't^r^et\vd /In-  squares . 

In  Figures  19  hnd  20  it  can  be  seen  that  thej|/e  jjases  range  from 
t;h3  equivalent  of  almost  perfect  .placement  to  almost  random 

‘ ■ . ...  -I  • \ . 

placement.  ' - 

Some  typical  fraction  damage  as  a function  of  distance 
examples  are  shown  in  Figure  22.  They  represent  a range  from 
very  steep  cutoffs  In  weapon  lethality  to  quite  gradual  ones. 

The  labels  give  an  indication  of  the  type  of  weapons  target  com- 
binations considered.  In  obtaining  these  curves  fraction  kill 
as  a function  of  distance  Is  obtained  by  combining  a function  of 
fraction  damage  as  a function  of  overpressure  with  a function 
of  overpressure  as  a function  of  distance.  The  fraction  of 
damage  as  a function  of  overpressure  is  assumed  to  be  a cumu- 
lative log  normal  curve  and  the  value  a in  Figure  22  is  the 
ratio  of  the  standard  deviation  of  the  cumulative  log  normal 
curve  to  the  mean  lethal  overpressure. 


Figure  22.  PROBABILITY  OF  KILL  AS  A FUNCTION  OF  NORMALIZED  DISTANCE  FOR 
SEVERAL  KILL  FUNCTIONS 
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